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We report a search for effects of large extra spatial dimensions in pp collisions at a center-of-mass 
energy of 1.8 TeV with the D0 detector, using events containing a pair of electrons or photons. The 
data are in good agreement with the expected background and do not exhibit evidence for large 
extra dimensions. We set the most restrictive lower limits to date, at the 95% confidence level, on 
the effective Planck scale between 1.0 TeV and 1.4 TeV for several formalisms and numbers of extra 
dimensions. 

PACS numbers: 04.50.+h, 04.80.Cc, 11.25.Mj, 13.85.Rm 



The possibility that the universe has more than three 
spatial dimensions has long been discussed [Q. Recent 
developments in string theory suggest that there could 
be up to seven additional spatial dimensions, compact- 
ified at very small distances, on the order of 10~^^ m. 
In a new model inspired by string theory, several of 
the compactified extra dimensions (ED) are suggested to 
be as large as 1 mm. These large ED are introduced 
to solve the hierarchy problem of the standard model 
(SM) by lowering the Planck scale, Mpi, to the TeV 
energy range. The ED compactification radius, R, de- 
pends on the number of extra dimensions (n) and on the 
effective Planck scale, Ms: R oc ;^(Mpi/Ms)2/" §. 
Since Newton's law of gravity would be modified in the 
presence of compactified extra dimensions for interaction 
distances below the size of the largest extra dimension, 
current gravitational observations rule out the case of a 
single large extra dimension. Recent preliminary results 
from gravity experiments at submillimeter distances 
as well as cosmological constraints from supernova cool- 
ing and distortion of cosmic diffuse gamma-radiation , 
indicate that the case of n = 2 is likely ruled out as well. 
However, for n > 3, the size of the ED becomes micro- 
scopic and therefore eludes the reach of direct gravita- 
tional measurements. Cosmological constraints are also 
weak in this case. Therefore, high energy colliders, ca- 
pable of probing very short distances, are crucial to test 
theories of large ED. In these theories, the effects of grav- 
ity are enhanced at high energies due to the accessibility 
of numerous excited states of the graviton (referred to 
as a Kaluza-Klein graviton, Gkk), corresponding to 
multiple winding modes of the graviton field around the 
compactified dimensions. Since gravitons couple to the 
energy-momentum tensor, they can be produced in any 
SM process. 

Large ED phenomenology at colliders has been studied 
m some detail § | |. One of the primary observ- 
able effects would be an apparent non-conservation of 
momentum caused by the direct emission of gravitons, 
which leave the three flat spatial dimensions. A typi- 
cal signature would be the production of a single jet or 
vector boson at large transverse momentum. The other 
observable effect would be anomalous difermion or dibo- 
son production at large invariant masses (M) via virtual 
graviton exchange. Direct graviton emission is expected 
to be suppressed by a factor (M/M5)"+^, while virtual 
graviton effects depend only weakly on the number of ex- 



tra dimensions ^. Virtual graviton production there- 
fore offers a potentially more sensitive way to search for 
manifestations of ED ||l^ . 

Both of the above effects have been sought at LEP , 
with lower limits on the effective Planck scale set on the 
order of 1 TeV. Virtual Gkk exchange was also sought 
at HERA 1 12 , but with less stringent limits. In this 
Letter, we report the results of the first specific search 
for ED at a hadron collider. The data are from the 
D0 experiment ||l^ at the Fermilab Tevatron, using 
proton-antiproton collisions at a center-of-mass energy 
of 1.8 TeV, and final states containing pairs of electrons 
or photons. We analyze the differential distribution |lj] 
of diclectron or diphoton events in terms of their invari- 
ant mass and the scattering angle in their center-of-mass 
frame {9*). These two variables completely define the 
2 — > 2 scattering processes, and their properties provide 
improved separation between the SM contributions and 
the effects of virtual graviton exchange when compared 
to analyses of invariant-mass distributions alone ||lj| . 

Several modifications to the method of Ref. ||l^ are in- 
troduced to optimize the sensitivity of the search. First, 
because the efficiency and resolution for high-energy elec- 
tromagnetic (EM) objects at D0 are superior to those for 
muons, we use only the diclectron and diphoton channels, 
with the sensitivity to ED coming primarily from the 
diphoton events iQ. Second, to optimize the efficiency 
for diphoton and diclectron selection, we eliminate track- 
ing requirements from electron and photon identification, 
thereby effectively combining them. Ignoring tracking 
information maximizes both the diclectron and diphoton 
efficiencies since neither electrons with unreconstructed 
tracks nor photons with matching tracks from conversion 
or random overlap are lost. In what follows, we refer to 
electron or photon objects that do not use tracking infor- 
mation as EM objects. The dominant background at high 
mass, where ED effects are enhanced, is the irreducible 
SM background from direct diphoton production, rather 
than instrumental background from misidentification of 
jets as EM objects (see, e.g., Ref. [|l^). Increasing the 
di-EM identification (ID) efficiency by using looser EM 
ID criteria therefore results in better sensitivity to ED. 

The search is based on the entire sample of data col- 
lected during 1992-1996 using a trigger that requires the 
presence of two EM objects and corresponds to an inte- 
grated luminosity of 126.8 ± 5.6 pb~^. 

Following offline reconstruction, we require each EM 
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object to: (i) deposit more than 95% of its energy in 
the EM calorimeter; (ii) have an energy deposition pat- 
tern consistent with that expected for an EM object; and 
(iii) be isolated ||l6|. The efficiency for these selections, 
based on Z ^ ee events, is 87 ± 2% per EM object. 
We also require (i) two EM objects with Et > 45 GeV 
in good fiducial regions of the detector {\r]d\ < 1.1 or 
1.5 < |?7d| < 2.5, where rjd is the pseudorapidity relative 
to the center of the detector |l^); (ii) missing transverse 
energy < 25 GeV; and (iii) no additional EM objects 
with Et > 5 GeV. The trigger is fully efhcient for this set 
of offline selections. The above selections, with the over- 
all efficiency of 79 ± 2% per event, define our base sam- 
ple, which contains 1282 candidate events. The error on 
the efficiency includes uncertainties due to the statistics 
and background parameterization in the Z — > ee sample. 
Other checks on possible variation of selection efficiency 
with energy, pseudorapidity, and invariant mass indicate 
that the efficiency is constant within uncertainties. 

To determine the hard-scattering vertex, we calcu- 
late the most probable direction of each EM object us- 
ing the transverse and longitudinal segmentation of the 
EM calorimeter Among all the reconstructed ver- 

tices in the event, we choose the one that best matches 
these directions. We choose this EM-object-based vertex- 
finding technique since it treats the diphoton and dielec- 
tron events in the same way. 

We model the effects of ED via the parton-level 
leading-order (LO) Monte Carlo (MC) generator of 
Ref. jl^ , augmented with a parametric simulation of the 
D0 detector. The simulation takes into account detector 
acceptance and resolution for the EM objects, smear- 
ing and misidentification of the primary interaction ver- 
tex, initial state radiation, and the effect of different par- 
ton distributions. We used leading order CTEQ4L0 Q 
distributions to estimate the nominal prediction. The 
parameters of the detector model are tuned using inde- 
pendent samples of collider data. The MC includes SM 
contributions {Z/j* and direct diphoton production), 
Kaluza-Klein graviton exchange diagrams, and their in- 
terference in di-EM object production. 

Since the parton-level generator involves only the 2 — > 
2 hard-scattering process, we model next-to-leading or- 
der (NLO) effects by adding a transverse momentum to 
the di-EM system, based on the transverse momentum 
spectrum of di-EM events observed in the data. In the 
presence of the NLO corrections, the scattering angle 9* 
is defined in the di-EM helicity frame, i.e., relative to the 
direction of the boost of the di-EM system. Since the 
parton-level cross section is calculated at LO, we account 
for NLO effects in the SM background by scaling the cross 
sections by a constant i^T- factor of 1.3 ||l^]. We assign a 
10% systematic uncertainty on the value of the if-factor 
to account for a possible growth of the if -factor at high 
mass. Because NLO corrections to the Kaluza-Klein di- 
agrams have not yet been calculated, we use the same 



constant if-factor for the signal. The it'-factor for gravi- 
ton exchange is expected to grow with invariant mass, 
similar to that for Z/^* exchange; consequently, our as- 
sumption tends to underestimate the ED contribution at 
high invariant mass and is conservative. 

The main SM sources of di-EM production. Droll- Yan 
and direct diphoton production, are included in the MC. 
Other SM sources, such as W'f, Zj, WW, Z — > rr, and 
ti production are negligible, as most have small cross 
section compared to that for Drell-Yan and direct dipho- 
ton processes, and are reduced further by the require- 
ments of low ]^rp and exactly two EM objects in the 
event. The only non-negligible sources of background 
arise from single-photon and dijet events, in which jet(s) 
are misidentified as EM objects. We estimate this back- 
ground using a separate data sample collected with a 
trigger requiring a single EM object. We require that a 
combination of a jet and the EM object satisfy all the 
requirements for signal, except that the EM ID require- 
ments are applied only to the EM object. We obtain 
the instrumental background by scaling the number of 
jet-EM combinations passing our selection criteria by the 
probability of a jet to be mistaken as an EM object, mea- 
sured to be 0.18 ± 0.04%, and, within the uncertainties, 
is independent of Et or rj |20 . The other source of in- 



strumental background, W-\- jets production with a jet 
misidentified as an EM object, is negligible due to the re- 
quirement of low $!rp. The total instrumental background 
in the base sample is 87 ± 22 events, which is less than 
7% of the dominant SM background. The sum of the SM 
background and the instrumental background reproduces 
the main kinematic characteristics of the base sample, as 
illustrated in Fig. ^ 

The cross section in the presence of large ED is given 
by i I, i: 



a 



dMdcos 6* 



fsM + fintVG + /kK?7g> 



(1) 



where /sm, /int, and /kk are functions of (M, cos^*) 
and denote the SM, interference, and Gkk terms. The 
effects of ED are parametrized via a single variable rjQ — 
J-/Mg, where ^ is a dimensionless parameter of order 
unity, reflecting the dependence of virtual Gkk exchange 
on the number of extra dimensions. Different formalisms 
use different definitions for T: 



T = 
T = 



1, (GRW [6]); 
/ log (S) ' " 



2 

n-2 ' 



2 

n > 2 



2A 2 

T = — = ±-, (Hcwett 

TT TT 



(2) 

(HLZ §); (3) 
)■ (4) 



Note that only within the HLZ formalism T depends 
explicitly on n. 

Figure |^ shows the two-dimensional distribution in M 
and |cos6'*| for data in (a), total background in (b). 
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-2 2 

Pseudorapidity, leading 



-2 2 

Pseudorapidity, second leading 




50 100 150 
E^, leading, GeV 



50 100 150 200 
E^, second leading, GeV 



FIG. 1: Comparison of data (points with error bars) and 
background predictions (solid histogram) for the pseudora- 
pidity and transverse energy of the two EM objects in the 
event. The shaded histogram shows the contribution from in- 
strumental background. The dips in pseudorapidity (defined 
relative to the primary interaction vertex) reflect the accep- 
tance criteria. 



and the sum of the background and the ED signal for 
i]c — ^ TeV^''^ in (c). The data agree well with the back- 
ground prediction and do not exhibit evidence for large 
ED, which would produce an excess of events at large 
mass and small values of |cos0*| (see Fig. ^). A com- 
parison of Figs. ||d and |^ shows that use of both M and 
cos 0* can improve the sensitivity over just one variable. 
The projection on cos 9* for high mass values in Fig. ^ 
shows the relative increase of signal over SM processes 
at small [cos 6**1. The use of both variables allows about 
10% improvement in the sensitivity to rjc- 

In the absence of evidence for ED, we set limits on 
the effective Planck scale. We perform a Bayesian fit of 
the sum of the cross section given by Eq. (|l|) and the 
instrumental background, to the data in the entire (M, 
cos 0*1) plane shown in Fig. ^, with r]G as free param- 
eter, with an assumed uniform prior distribution. The 
systematic uncertainties on signal and background in the 
fit include systematics of the if -factor (10%), choice of 
parton distribution functions (5%), integrated luminos- 
ity (4%), EM ID efficiency (3%), and the uncertainty on 
the instrumental background (25%). 

The best estimate of the parameter rjc is consistent 
with the SM value of 77^ = 0, and the one-sided 95% CL 
limits on r]G are: 



riG < 0.46 TeV"'' (?7g > 0) 
r]G > -0.60 TeV"* {tjg < 0), 



(5) 
(6) 




SM+ED signal, r|(, = 1 TeV 



0.25 0.5 0.75 1 
|cos(e*)|, M > 250 GeV 



FIG. 2: Two-dimensional distributions in di-EM mass and 
j cos 6**1 for: (a) data, (b) background, (c) background and 
ED signal for 770 = 1 TeV~^, and (d) | cos6*| distribution for 
events with M > 250 GeV, where the filled circles correspond 
to the data, instrumental background is shown shaded, the 
entire background from SM sources is given by the solid line, 
and the dotted line corresponds to the sum of SM and ED for 
VG = 1 TeV-*. 



in good agreement with the expected sensitivity to tjg, 
as obtained in an ensemble of MC trials (0.44 TeV^^ for 
Vg > 0). We can express these results in terms of lim- 
its on the effective Planck scale for the three formalisms 
of Eqs. (|)-(|). In the formalism of Ref. §, both signs 
of rye are possible and therefore both limits (^) and (||) 
are relevant. In the other two formalisms, rjG is always 
positive, and only the first limit is relevant. For the HLZ 
formalism, the case of n = 2 is special since and there- 
fore rjG, depend on M. To relate rjG to AIs for n — 2, 
we use an average M for the Gkk term at the Tevatron 
of (0.6 TeV)^ [Q. Limits for different formalisms and 
for different numbers of extra dimensions are given in 
Table |. They correspond to the ED compactification 
radius ranging from R < 0.3 mm (n = 2) to i? < 2 fm 
(n = 7). 

TABLE I: Lower limits at 95% CL on the efi'ective Planck 
scale, Ms, in TeV. 



GRW 


6 


HLZ 






Hewett 1 




1.2 


n—2 n=3 n=4 n=5 n—6 n—7 


A = +1 A = -1 


1.4 1.4 1.2 1.1 1.0 1.0 


1.1 1.0 



In summary, we have performed the first search for 
large extra spatial dimensions at hadron colliders by look- 
ing for effects of virtual Kaluza-Klein gravitons in the 
production of dielectrons and diphotons at high energies. 
No evidence is found for large ED, and we set a model- 
independent 95% CL upper limit of 0.46 TeV""' on the 
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parameter r]a that describes the size of the ED contri- 
bution. This result corresponds to Hmits on the effective 
Plai iclt scale ranging between 1.0 and l.d ToV for several 



forn ialisins and numbers of ED. Those are the most re 



strictive limits on large ED to date, and are complemen- 
tary to analogous limits from LEP that probe a different 
range of invariant masses. 
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